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ABSTRACT 

The Local Group contains some dwarf galaxies that apparently formed all their stars early, lost 
their gas to winds or ionization, and survive today as "fossils" from the epoch of reionization. This 
study presents new models of these objects based on the hierarchical chemical evolution framework 
of Tumlinson (2006). The model accurately reproduces the observed luminosity-metallicity relation 
of fossil dwarf galaxies with minimal free parameters. When calibrated to this relation, the models 
show that small dark matter "minihalos" formed 2 — 8% of their baryonic mass into stars prior during 
reionization. By tracking the chemical enrichment of these early halos the models specify that metal- 
free first stars contributed approximately 5 — 10% of the ionizing photons generated by these small 
halos and so did not dominate reionization. Models that allow for larger relative contributions from 
metal-free stars may not generate enough total ionizing photons for early reionization. As significant 
star formation in early minihalos is a requirement of many successful models for the IGM reionization, 
these models can be considered to pass a key test of their validity. It appears we have as much to 
learn about reionization from the universe at z = as at z = 6. 

Subject headings: galaxiesiformation - galaxies: evolution - Galaxy: format ion - Galaxy:evolution — 
stars:abundances - stars:mass function - cosmology:theory 



1. INTRODUCTION 

The reionization of the intergalactic medium was a ma- 
jor change in the physical state of most of the baryons 
in the Universe, potentially a major influence on the for- 
mation and development of galaxies, and may be inti- 
mately connected with the first generations of metal- 
free stars. Observations by the Wilkinson Microwave 
Anisotropy Probe ( WMAP) of the integrated Thomson 
scattering optical depth to the cosmic microwave back- 
ground, Tes = O.ITIq'qy, and the Gunn-Peterson ab- 
sorption trough in spectra of high-z QSOs found by the 
Sloan Digital Sky Survey (SDSS), together suggest that 
reionization had a complicated history in the interval 
z = 6 — 20, which spans < 1 Gyr of cosmic timc"'^. 

Theoretical studies have attempted to explain these 
data in terms of semi-analytic reionization models that 
calculate the collapse of dark matter halos and track the 
growth of cosmological H II regions (Venkatesan, Tum- 
hnson, & ShuU 2003; Haiman & Holder 2003; Cen 2003; 
Tumhnson, Venkatesan, & ShuU 2004, hereafter TVS04). 
These studies vary in their particulars, but they gener- 
ally include assumptions about or ab initio calculations 
of three important physical quantities that have so far 
been unconstrained by observation. The first of these is 
/* , the efficiency of converting virialized baryons within 
dark-matter halos into stars. This factor may depend on 
such local conditions as gas density, metallicity, and radi- 
ation feedback in unknown fashion, and as /* is difficult 
to calculate a priori, it is generally assumed. The sec- 
ond uncertain parameter is 70, the time-integrated ion- 
izing photon production per baryon in stars. This quan- 
tity is essentially completely determined by the initial 
mass function (IMF) at a given metallicity, and ranges 

^ I assume the cosmological parameters derived from first-year 
WMAP data by Spergel et al. (2003); = 0.27, Ha = 0.73, 

Ut = 0.044, and Ho = 71 km s"! Mpc"!. 



from 7o = 4000 for a Salpeter IMF at solar metallicity to 
^ 50000—100000 for an IMF dominated by massive stars, 
as expected at low metallicity. TVS04 used detailed 
chemical abundances in metal-poor Galactic halo stars 
have favored massive star IMFs with 7 ~ 50000 — 80000 
for metal-free stars. Finally, the escape fraction of ioniz- 
ing photons from their parent galaxies is denoted by /g^^ 
and has been estimated at 0.1 for bright galaxies in the 
local universe, but is unconstrained for small galaxies in 
the high-redshift universe. Thus, the three critical model 
inputs for reionization are all uncertain. 

Theoretical studies attempting to explain the reion- 
ization data have varied in their specific treatment of 
these uncertain parameters, specifically with regard to 
their assumptions about "minihalos" . TVS04 achieved 
Tes = 0.12 with = 0.05, fesc = 0.05, and an IMF with 
10 - 100 Mq, but did not vary /* and fesc- Haiman & 
Holder (2003) combined the three factors into a single 
efficiency parameter, Chaio = /*7o/esc, and varied this to 
achieve reasonable fits to the WMAP and Gunn-Peterson 
data. Cen (2003) studied the feedback effects of X-ray 
emission from early SN remnants and black hole accre- 
tion and found that low-mass halos could form stars ef- 
ficiently enough to reionize the universe completely by 
z = 17. Though they differ in their details, the theoreti- 
cal studies cited above generally agree that dark matter 
"minihalos" (with virial temperatures T = 10'^ — 10^ K, 
or M = 10^ - 10*^ Mq at z = 6 - 15), which could prob- 
ably cool only by molecular hydrogen (H2), are essential 
to explaining the WMAP data on Tes- The reasons are 
the cumulative nature of Tes and the bottom-up nature 
of structure formation. If the universe is completely ion- 
ized only at z < 10, Tes ~ 0.08. High levels of ionization 
at z > 10 are therefore necessary to match the WMAP 
constraint, which corresponds to complete ionization to 
z — 17. In the range z = 10 — 17, minihalo hold most 
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of the collapsed baryons in the ACDM cosmology. At 
z — 15, 3.3% of all baryons reside in collapsed halos 
with Tmr > 10^ K, but only 0.7% are in objects with 
Tvir > 10^ K. If ^ 80% of the baryons do not partic- 
ipate in star formation, extreme assumptions about 70 
and /esc would be necessary for = 0.17 or the models 
may fail altogether. 

Dark-matter "minihalos" apparently need to form 
stars with efficiency /* > 0.01 to match the key data on 
reionization. However, these objects have the most un- 
certain physics of gas cooling and star formation and are 
likely to be too faint to be constrained by independent 
observation at high redshift. Some hydrodynamical sim- 
ulations (Abel, Bryan, & Norman 2000, 2002; Bromm, 
Coppi, & Larson 1999, 2002) have followed H2 cooling in 
these objects and determined that they likely form only 
very massive stars (VMSs; M ^ 100 Mq) with low ef- 
ficiency (/* ~ 0.001 — 0.01) and then rapidly dissociate 
the H2 in their neighboring halos, thereby suppressing 
further the average minihalo efficiency. Other authors 
(Ricotti et al. 2002; Cen 2003) have argued that the 
feedback of soft-UV and X-ray photons is likely to be 
positive, such that minihalos can experience significant 
star formation over long periods. In light of all these re- 
sults, efficient star formation in these halos can be taken 
as a prediction of successful reionization models, most of 
which would fail without the participation of minihalos. 

This study extends previous work (particularly by Ri- 
cotti & Gnedin 2005, hereafter RG05) to argue that these 
predictions can be tested already using dwarf galaxies 
in the Local Group (LG), and that these dwarf galax- 
ies likely formed as "minihalos" at high redshift and 
contributed significantly to the reionization of the IGM. 
Section 2 of this Letter briefly describes the framework 
of Tumlinson (2006), concentrating on the modifications 
needed to study LG dwarfs. Section 3 reviews the data on 
Local Group dwarfs and presents models for their lumi- 
nosity, metallicity, and star formation efficiency. Section 
4 summarizes the conclusions. 

2. METHODS 

Tumlinson (2006, hereafter Paper I) presented a new 
chemical evolution code that works both within the hier- 
archical context of galaxy formation and in the stochastic 
limit of low-metallicity Galactic evolution. Paper I fo- 
cused on primordial IMF constraints using Galactic halo 
stars, but the framework is adaptable to model halos at 
arbitrary mass and redshift and so is well suited to mod- 
eling pre-Galactic halos that form before reionization. 

The framework uses the common technique of halo 
merger trees (Somerville & Kollatt 1999) to decompose 
the Galaxy into its precursor halos working backward in 
time. It then calculates the history of star formation and 
chemical enrichment in these objects working forward 
in time, keeping track of all individual metal-producing 
supernovae and assigning metallicity to new star for- 
mation stochastically from all prior generations. Stars 
are formed with a constant efficiency per timestep, e,, 
such that the mass formed into stars = e^Mgas^t in 
time interval At. At metallicities below the critical value 
for normal star formation, Zcrit = 1O~'^'^.Z0 (Schneider 
et al. 2002; Santoro & ShuU 2006), the IMF is a log- 
normal function with characteristic mass nic and width 
(TjMF- Above Zcrit ^ a power-law IMF with Salpeter slope 



{a = -2.35) is assumed, with 0.5 - 140 Mq. IMFs with 
higher rric or steeper a are more "top-heavy" , while Zcrit 
effectively sets the time over which the top-heavy IMF 
persists. In the small halos of interest, metallicity is cor- 
related with time but does not increase monotonically 
over short times. 

The Paper I method is changed slightly to model dwarf 
galaxies in the Local Group. The model dwarf galaxies 
are like their real-world counterparts in that they virial- 
ized early, formed stars for a short period, and then lost 
their gas, were ionized by internal or external sources, 
or otherwise had their star formation truncated at some 
redshift, Zend, the redshift of the root halo in the tree. 
After this redshift they evolve no further. In the real 
universe, these dwarf galaxies may continue to evolve 
dynamically by accretion into larger objects, tidal strip- 
ping, or they may remain in isolation, but they do not 
appear to evolve chemically beyond the loss of their gas. 
A large grid of models was produced with Mh = 10^ — 10^ 
Mq, e, = 0.1 - 3.0 X 10-1° yr-i, and Ze„d = 6-16. The 
output quantities of interest are sensitive to these quan- 
tities but not to the choice of IMF below Zcrit or to the 
interstellar mixing rates. 

3. LOCAL GROUP DWARFS AND REIONIZATION 

The roughly 40 known Local Group dwarf galaxies 
span a range of types (from spheroidal to elliptical to 
irregular), have generally low metallicity, and dynamical 
mass from 10° - 10^ Mq (Mateo 1998). Prompted by the 
old age of many dSph populations, and by the discrep- 
ancy between CDM predictions for dwarf galaxy num- 
bers and the known population, recent theoretical stud- 
ies have posited a close connection between reionization 
and the early evolution of LG dwarf galaxies. Bullock, 
Kravtsov, & Weinberg (2001) suggested that the suppres- 
sion of star formation in small halos by ionizing photons 
at z > 6 could explain the relatively small number of LG 
dwarfs compared with the CDM prediction. 

RG05 examined this "reionization hypothesis" in hy- 
drodynamical simulations that included II2 cooling and 
radiative feedback. They were able to reproduce many 
of the observed properties of the LG dwarfs in simula- 
tions that stopped at z = 8.3, with no further evolution 
in the galaxies of interest. They proposed a classifica- 
tion scheme in which LG dwarfs are divided into three 
categories: "survivors" massive enough to retain their 
gas and continue forming stars robustly after reioniza- 
tion, "polluted fossils" that were affected by reionization 
but managed to continue forming stars at a modest level, 
and "true fossils" , which formed most of their stars be- 
fore reionization and essentially none afterward. I adopt 
this scheme here and use the "true fossils" , those galax- 
ies with exclusively old stellar populations, roughly 1% 
solar metallicity, and My > —12, to constrain the star 
formation efhciency of small halos at high redshift and 
assess the star formation efficiency of minihalos. 

The first result of this study is that the observed 
luminosity-metallicity (L — Z) relation of LG dwarfs is 
reproduced exactly with minimal tuning of the mod- 
els. This relation appears in Figure ^ where a selected 
group of LG dwarfs is plotted in the style of Caldwell 
(1999). Survivors and polluted fossils are marked in gray 
type; true fossils are italicized. Mean metallicities are 
taken from Mateo (1998), except for And V (Zucker et 
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al. 2004) and Ursa Major (Willman et al. 2005). Ab- 
solute V-band magnitudes for z = are obtained by 
interpolation of the tabulated isochrones for 1/50 and 
1/200 solar metallicity from Girardi et al. (2002) and 
applied individually to stars that reside in the final halo 
at Zend- The points mark the theoretical calculations for 
Mh = 10^10^5 X 10^2 X 10^, and 10^ Mq, from left 
to right; the models are not intended to perfectly match 
the survivor and polluted fossils that have experiences 
recent star formation (in gray). Within each constant- 
mass sequence, the time-averaged star formation param- 
eter e, = 1 — 3 X 10^^° and Zend = 6—10. An increase 
in e* or decrease in Ze_nd moves points up and to the left 
in the constant-mass sequences. For the same e^, and 
Zend combination, the variation in mass traces out a line 
of the same slope as the observed trend. We can there- 
fore interpret the observed L — Z relation as a variation 
in mass with scatter introduced by small variations in 
the star formation histories of individual objects. This 
concordance was achieved with minimal tuning; it was 
necessary to vary only with the same mean and in 
the same range found in Paper I to match the Galactic 
halo metallicity distribution. I conclude that the sim- 
ple description of LG dwarf satellites as fossils from be- 
fore reionization is accurate and that a more complicated 
model is unnecessary. 

The second result of this study is the strong correla- 
tion between total star formation efficiency, /* , and mean 
metallicity of the resulting dwarf galaxy for thousands 
of model halos (Figure . Because metallicity and ef- 
ficiency are effectively mass-normalized quantities, this 
correlation is insensitive to halo mass. Only at the low- 
est mass plotted here {Mh = 10^ Mq) does the relation 
depart to slightly higher efficiency for a given metallicity, 
with increased scatter. This relation can now be used to 
determine the star formation efficiency of Local Group 
"fossil" dwarfs. 

Estimates for f^, in the LG fossil dwarfs are obtained 
by varying the assumed total halo mass, Mh and time 
averaged star formation parameter, e,, to obtain best 
fits to the observed My and {[Fe/H]), assuming the er- 
rors in these quantities tabulated by Mateo (1998). The 
optimization is done with a Monte Carlo Markov chain 
maximum likelihood estimator so that confidence inter- 
vals are also obtained. Figure 13 shows the results for 
total star formation efhciency /* compared with Tyir de- 
rived from fits to the observed L-Z relation. For the 
model "true fossils" , /* ranges from 2 — 8%, while the 
best-fitting halo masses clearly indicate that these galax- 
ies lie in the regime of cosmological "minihalos" at high 
redshift. There is also tentative evidence for a positive 
correlation of /* with T^r, at 90% significance; a lin- 
ear best fit is shown in the dotted line. These results 
assume a fixed z^nd = 6, but there are no substantial 
changes to the individual values or trend if z^^id — 10. 
The width of the model metallicity distribution function 
(MDF) does vary with Zend', shorter total star formation 
times (higher z^nd) yield narrower distributions. Halos 
that form stars until z = 6 have broader distributions 
than those that had their star formation truncated ear- 
lier, but the existing observations of a{[Fe/H]) are too 
imprecise to use this as a sensitive indicator of the dura- 
tion of their bursts. Constraints on models of the detailed 
star formation and feedback mechanisms in these objects 



would benefit from the more sensitive and precise obser- 
vations of MDF and relative chemical abundances that 
future large telescopes could provide. 

Based on these results I conclude that minihalos con- 
tribute significantly to the reionization of the IGM at 
z > 6. However, we cannot be sure that the minihalos 
we see today as fossil dwarf galaxies in the LG are typical 
of minihalos that formed in this region of space, or that 
the LG population is representative of the minihalo pop- 
ulation as a whole. Such a concern could be answered by 
detailed examination of dwarf galaxies outside the Lo- 
cal Group, but this task is probably too challenging for 
observations in the near future. Future theoretical work 
should attempt a full chemo-dynamical treatment of the 
Local Group to help assess the bias of the LG dwarfs in 
mass or metallicity, which may change their place in the 
scheme of reionization. 

Total masses inferred from the L-Z relation in the pro- 
cess of constructing Figure 01 are compared to their dy- 
namically estimated masses in Figure 01 where the best- 
fit Mh is plotted against the measured Mtot for the nine 
fossils with rotation-curve mass estimates (Mateo 1998). 
The models generally best match the observed L-Z re- 
lation for masses similar to the dynamical estimates (a 
1:1 ratio is marked with the dashed line). Thus it does 
not seem necessary to invoke significant mass loss from 
the "true fossils" to explain any of their observed prop- 
erties, except in a few cases. By this measure. Sculptor, 
Antlia, and And VI appear to have lost significant mass 
in this interval. This apparent mass loss, compared with 
other galaxies, can be taken as a prediction of the mod- 
els presented here that may be testable by more detailed 
examination of their stellar populations. 

The merger-tree models keep individual account of all 
supernova-forming stars, so the distribution of stellar 
metallicities is known exactly at all times. The models 
specify the fraction of ionizing photons produced by stars 
of different metallicities. Figure shows the fraction 
of ionizing photons from metal- free stars, Fz=o, plotted 
against the total ionizing photon budget per baryon in 
the final halo {N.y w /»7o > 1000 for the t^s ^ 0.1 models 
of TVSG4) for four of the five IMF test cases of Tumlin- 
son (2006). For z^nd = 6 (open circles), the mean Fz=o 
is « 5% for IMFs A and D. Results for IMF B are similar 
to case A and are not shown in the figure. Cases C and 
E are more top-heavy and give Fz=o ~ 0.5 with a long 
tail to high Fz=o- Paper I showed that cases A, B, and 
D, provided the best overall fit to the data on Galactic 
chemical evolution, while cases C and E generally over- 
produce very massive stars that are inconsistent with the 
abundance data. If z^nd = 15 (open squares), the relative 
contribution of metal-free stars is larger, Fz=o ~ 0.1—0.5 
for cases A and D, but here stars are formed for a shorter 
total time and so have reduced ^ 200. As reioniza- 
tion models typically require N~f > 1000 from minihalos 
(TVSG4, HHG3) to achieve t^s > 0.1, it seems that a suc- 
cessful model cannot reproduce the WMAP t^s and also 
have metal-free stars as the primarily agents of reioniza- 
tion. If the minihalos that survive are typical of those 
that formed during the epoch of reionization, it appears 
that metal-free stars started but could not have com- 
pleted this critical transition in the state of matter in 
the IGM. 
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4. CONCLUSIONS 

Based on halo-merger-tree chemical evolution models 
of Local Group "fossil" dwarfs, I conclude that: 

1. Models of small dark matter halos forming stars 
during reionization accurately reproduce the ob- 
served luminosity-metallicity relation for LG dwarf 
galaxies with the same model parameters that also 
match the Galactic stellar halo. 

2. The strong correlation between mean metallicity 

and total star formation efficiency in "true fossils" 
implies total star formation efficiency /* = 2 — 8% 
for halos with mass Mh = 10'' - 10^ Mq, or T^ir = 
2000 - 7000 K at z = 6 - 15. 

4. By tracking chemical enrichment inside these small 
halos, the models also specify that metal-free stars 

contribute approximately 5 — 10 % of the total ion- 
izing photon budget from minihalos to z = 6. Thus 
metal-free stars probably did not provide the ma- 
jority of ionizing photons to the IGM. 

I conclude that dwarf galaxies that form all their star 
early and survive into the present-day Local Group pro- 
vide sensitive indicators of the star formation history of 
low-mass "minihalos" during the epoch of reionization. 
By comparison with the results of RG05, which stim- 
ulated this work, my results more accurately describe 
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the present-day luminosity-metallicity relation and in- 
dependently derive the total star formation efficiency 
and first-stars contribution to the ionizing budget. Con- 
versely, these models lack the sophisticated treatment of 
radiative feedback and three-dimensional hydrodynami- 
cal evolution that simulations allow and do not account 
for possible biases in the surviving dwarf galaxy popu- 
lations. The two approaches therefore complement each 
other and show how Local Group fossils provide poten- 
tially critical tests of competing scenarios for reioniza- 
tion. These results suggest that no reionization model 
should be considered successful unless it correctly incor- 
porates the metallicity and age of old LG stellar popu- 
lations into the history of this critical epoch. Because 
LG dwarf galaxies provide important information that is 
not available directly at high redshift, we may have as 
much to learn about reionization from z = 0, where we 
can study it residues, as at z = 6, where it ended. 
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Fig. 1. — Luminosity-metallicity {L — Z) relation for model dwarf galaxies on constant-mass sequences with Mj, = 10^, 10*, 5 X 10^,2 X 
10^, 10^ Mq from left to right. The star formation parameter e ranges from 1 — 3 X 10"^" yr~^ from the bottom to top of the constant-mass 
sequences. The Ursa Major dwarf (Willman et al. 2005) has only a metallicity range. 
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{[Fe/H]) 

Fig. 2. — The total star formation efficiency, /«, in small halos of varying mass that cease forming stars at 2e„j = 6. The symbols are 
the same as in Figure The correlation of /* with mean metallicity is strong and insensitive to M/^ and Zend- 
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Fig. 3. — The total star formation efficiency /» for the LG true fossils, derived from maximum likelihood fits to the observed mean 
metallicity and luminosity. Tyir is evaluated at z = 6 using the best-fitting halo mass, Mf^, and assuming this is unchanged since z = 6. 
The error bars mark the 2a (95%) confidence intervals derived from Markov chain searches of (Af^, e*) parameter space. The dotted line 
marks a tentative linear fit to (logT^ir, f»). 
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Fig. 4. — Comparison of total halo mass, M^, inferred from halo merger tree models with dynamical masses, Mtot- Error bars show 2cr 
(95%) confidence intervals. 




Fig. 5. — Ionizing photon fraction from metal-free stars, -Fz=Oi versus total ionizing photon production, N-y, for four of the IMF test 
cases in Tumlinson (2006), e.g. IMF A has rric = 10 Mq and (Jimf = l-O- 



